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Abstract 

We present some applications of the renormalized RPA in bosonic field the- 
ories. We first present some developments for the explicit calculation of the 
total energy in <I> 4 theory and discuss its phase structure in 1 + 1 dimensions. 
We also demonstrate that the Goldstone theorem is satisfied in the O(N) model 
within the renormalized RPA. 



1 Introduction 

The application to quantum field theories of non-perturbative methods used in the 
nuclear many-body problem has recently given rise to numerous promising works. 
One central motivation is to obtain tools to describe matter made of strongly in- 
teracting hadrons in presence of broken symmetries such as chiral symmetry. In 
particular the well-known RPA method, which has been originally developed in the 
context of condensed matter physics, has been recently applied to study a bosonic 
0(N) model (i.e. the linear sigma model). It has been demonstrated [2] that the 
standard RPA is able to restore the Goldstone theorem which is violated at the level 
of the usual variational Gaussian approximation [3J HJ [S| . Although this result can 
be seen as a major success, the RPA method in its standard form possesses some 
weak points. In particular it has the tendency to overestimate the attractive correla- 
tion energy, at least in examples of nuclear physics. In a recent paper JBJ, hereafter 
referred as I, we have developed the formalism of a superior version of the RPA, 
namely the renormalized RPA (r-RPA) in the particular context of A$ 4 theory. As 
shown in I, one important merit of the r-RPA is to cure the instability problem 
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appearing in the standard RPA in 1 + 1 dimensions. The purpose of this paper is to 
present a more detailed study of the phase structure of this theory within the r-RPA 
approach. Besides we introduce improvements of the former calculation and obtain 
a second order phase transition for this model. We also present the r-RPA method 
for the O(N) model, demonstrating that the Goldstone theorem (massless pions) is 
also satisfied at this level. Since all the details of the approach have been given in 
I, we limit ourselves to the strict minimum for what concerns the formalism. 

2 The r-RPA in $ 4 theory 

We consider the Lagrangian density: 

C = ±ff>*{x,t)d lt *{x,t) - ±fi&{x,t) - ^\x,t) (1) 

where /Xq is a constant and the bare coupling constant b = A/6 is positive for reasons 
of stability. We decompose the scalar field <J>(x,t) in a classical part or condensate 
s and a fluctuating piece (p(x,t): 

$(x,t) = <f>(x,t) + s, s = ($(x,t)). (2) 

The presence of the condensate s indicates a spontaneous breaking of the underlying 
<3? — > — symmetry. Introducing the conjugate field H(x), one obtains for the 
Hamiltonian (in d + 1 dimensions): 

H= jd d x j^ 2 + A, 4 + (^ 8 + *{x) 

+ ^[n 2 (x) + {di4>f{x) + (nl +^s 2 )4> 2 (x)] 

+ T 03(x) + 2l^ (x) }- (3) 

Putting the system in a large box of volume V = L d , it is convenient to work in 
momentum space and to expand the fields according to: 

^) = ^E e ^>^) . n(*) = --^J>^iv(i), (4) 

q q 

i.e., in terms of creation and annihilation operators obeying the standard canonical 
commutation relations: 



The choice of the basis (i.e. the choice of the Kg) will come out as a part of the RPA 
solution. 



In I |H] we have explicitly solved the r-RPA problem, using the Green's function 
method, taking into account one-particle bt and two-particle b^b^,, b^b_^/ excitation 
operators. We refer the reader to I for the detailed derivation and we only quote 
here the main results (we also give in the appendix the results for the 1-1 1-2 and 
2-2 Green's functions). 

It is convenient to choose the basis = where e$ is the generalized mean-field 
energy, solution of the gap equation (\H,b^b_^\) = 0. This gap equation reads: 

2 -2 , 2,^2, ^ / i2\ -2 , 2 i n\ 

£q = Q +^o+2 s +2^)r = ( i + m ( 6 ) 

where ((f) 2 ) r = (^/V)Ylq(^ + ^\^q)/^ £ q ls * ne self-consistent scalar density. To 
obtain the standard RPA case, one simply has to replace the self-consistent scalar 
density (4> 2 )r by the Gaussian one (4> 2 ) e i- e - the expectation value is calculated on 
the vacuum of particles having the energy and such that (btb$) e = 0. In one 
spatial dimension, the generalized mean-field mass m is rendered finite by a simple 
mass renormalization: 

2 2 _L 6 2 _L 6 (lA?\ / +A * 1 ^ m 

m =" + r + 2 {<*>*- L sfv?T7j <7) 

where \x is the renormalized bare mass of the theory. 

The r-RPA single particle propagator has been derived in I and is given by: 



G(E,P) = ( E 2 - e\ - E(E,P)) with E(E,P) = ^- ~\ J/J ^ ■ (8) 



b 2 s 2 I{E,P) 

imp) 



The two-particle loop I(E, P) has the explicit form given by eq. 1)46 [I in the appendix. 
It explicitly depends on the "occupation number" N$= (4>q4>cr)R which constitutes 
the remaining problem to solve. One serious difficulty is that covariance is lost in the 
r-RPA in the sense that the loop integral I(E, P) and consequently the mass operator 
T,(E, P) depends separately on E and P due to the presence of the density N in its 
expression. This is certainly a weakness of the present approach (see discussion in 
I). One natural possibility to recover covariance consists in imposing that the correct 
I(E, P) is obtained through its center of mass (CM) expression according to: 

I(E,P) = I(E 2 - P 2 )= f j^-j . (9) 

1 ' ' K ' J (27T) d £2 _ p2 _ 4£ 2 +irj ^ 



The densities Mt can be calculated self-consistently using the spectral theorem: 



Mp = J l ^e^ + G(E,P). (10) 

In the quasi-particle approximation used in I, the solution is A/p = l/2ftp where 
ftp = VM 2 + P 2 is the energy of the one-particle RPA mode. As explained in 
details in I, the problem reduces to determine the pole of the one-particle propagator. 
One important result of I is linked to the fact that the instability present in 1 + 1 
dimensions at the level of the standard RPA (imaginary solution for the mass M of 
the RPA mode) just disappears in r-RPA. 

3 r-RPA correlation energy in <J> 4 theory 

The energy of the system is calculated from the various Green's functions using 
the spectral theorem. However, as is well-known, this cannot be done directly from 
the RPA results since important correlations would be lacking at the level of the 
expectation value of the kinetic energy. To solve this problem we have generalized 
the so called "charging formula" method to the r-RPA case. We decompose the 
Hamiltonian in two pieces Hq and Hint- 

H = V Q/xg s 2 + A /) + \ (iT U\ + 0\ <Ml) + H 3 + H 4 = H + H int (11) 

where #3 and H4 are the 3- and 4-particle pieces of the Hamiltonian (last line of 
eq. ©)• Hq has a form of a free Hamiltonian for quasi-particles with the generalized 
mean- field mass m (eq. (J7J)): 

Ho = £ + ^ (: nx nt i -e +4- 4>1 4 -e) ■ (12) 

where Eq is the generalized mean-field vacuum energy given in eq. (61) of I. As 
explained in subsection 4.2 of I, the interacting Hamiltonian is thus: 

H mt = H z + H A - h - {<j?) R £ : 0a </>{ : E - V ^ (<p 2 ) 2 . (13) 

We also introduce an auxiliary Hamiltonian: 

H'(p) = H + P H int , H'(p = 1) = H . (14) 

The price to pay is to solve the r-RPA problem (in practice the calculation of the 
commutators and double-commutators entering the RPA equations) for the H'(p) 



Hamiltonian. As explained in I, this can be done by making the following modifica- 
tions in the corresponding r-RPA problem for H: 



H - H\p) 



b 



+ ^ ((<P 2 )rp - (<P 2 )r) 



H 3 + H± p(H 3 + H 4 ) (15) 

where {4> 2 )r p is the self-consistent scalar density in the correlated RPA ground state 
of H'(p). We now employ the charging formula to calculate the ground-state energy 
as a function of the condensate s, i.e., the effective potential needed to study the 
phase structure of the theory: 



E RPA = E +[ ^( P H int ) p . (16) 
Jo P 

Using Wick theorem with respect to the vacuum of the quasi-particle with energies 
s p , the correlated part can be rewritten as: 

( P H mt ) p = (pH 3 ) p + (p:H 4 : ep ) p -V?±((<l?) ep -(<!?) e ) 2 

-V ?± ({ <?)r -{<?),,)({ 4?)* • (17) 

In this formula (4> 2 }r is as before the self-consistent scalar density of the original H 
whereas ( 4> 2 )r p corresponds to the same quantity in the H'(p) problem. {4> 2 ) £ is the 
scalar density in the generalized mean field vacuum (vacuum of quasi-particles with 
energy e^) in the H problem and {4> 2 ) £p corresponds to the equivalent quantity for the 
H'(p) Hamiltonian. The expectation values (pH 3 ) p and (p : H4 : ) p are calculated 
by using the spectral theorem applied to the 1-2 and 2-2 Green's functions relative to 
the r-RPA H'{p) problem. These are actually the main contributions noted. cow 
and corr of the correlation energy. However, at this level, we would like to precise 
one point. In the previous article I, we have taken into account the term: 



rr = 1 *£ 

p 



v p 4{{ 4>X P -{<P 2 )s) 2 



-V^(( <P 2 )R-(<P 2 )s P ) (( t 2 ) Rp -(<P 2 )s) 



(18) 



(called factorized term thereafter) in the following manner. The self-consistent scalar 
density in the p problem was calculated by using the spectral theorem: 



4 X 



In this article, we calculate the difference of scalar densities for the various masses 
with the following formula (A — ► +00): 



<A0 2 ) 
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We will show in the section about numerical results that this method gives better 
results. In particular, with the non covariant terms (see discussion below) we obtain 
this very important new result : the Simmon-Griffith theorem (which states that 
the order of the transition in the $ 4 model cannot be of the first order) is satisfied. 

Two problems arise in the calculation of the effective potential. First, there is 
an ambiguity in the calculation of the expectation value (pH^) p which was not 
addressed in I and second, non covariant terms appear, even when the covariance 
is forced in the loop integral I p (E, P). We will show below that the first ambiguity 
gives a numerically negligible effect and concentrate mainly on the non covariance 
problem. We notice that we have combinations of Green's function (see eqs. (|52|) and 
1|53|) in the appendix) which are explicitly covariant if the two-particle loop integral 
is covariant. Hence, we can decompose the 3-particle and 4-particle correlation 
energies in covariant pieces (which correspond to the expression (|52|) and (|53|)) and 
non covariant pieces (the remaining terms) according to: 



zp _ tp(3c) , K (3nc) , 7T.(4c) , pUnc) 

The 4-body pieces have been given in I with the result : 



E, 



(4c) 



E, 
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(4nc) 



V 



1 dp 

P 

1 dp 



dP 



idE 



(2ir) d J (2tt) 



e lE ^ + I 2 JE,P)F p (E,P), 



(20) 
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with FJE,P) 
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f idE 


(2vr) d . 
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p 2 b 2 



e iE v + 4 2 ^ ^ p p ^ ; (22) 



p 2 b 3 s 3 G P {E, P) 
l-fl p (E,P) U-^T^pf 



+ 



(23) 



where the indices p mean that the quantities are related to the H'(p) problem. The 
non covariance comes from the presence of the loop integral: 



I^(E,P) 



dk 



(27T) 



2 e\ e 2 



E — £\ — £2 + irj 



(24) 



which is not covariant (i.e., it depends separately on E and P) even if covariance 
is forced by taking as before the CM expression. We have neglected E^trr^ (which 
vanishes to leading order in the interaction) in our previous numerical estimate 
in I. However this contribution, although relatively small, turns out to be very 
important for the precise nature of the phase transition. In other words i?( 4nc ) is 
very important to reproduce good numerical results. Lets us come to the calculation 
of the expectation value of #3: 

(H 3 ) = - lL— 5 1+2+3 (((b\ b\ + 6_! 6_ 2 ) (6_ 3 + 4) 



+ 2(6146-3 + b\ 6-26-3)) )• (25) 



The first line generates the covariant contribution: 

E%£ [Up f dP f idE lErj+ p 2 b 2 s 2 G p {E,P)I p (E,P) 



V Jo P J (27r) d J (2tt) 6 l-e+I p (E,P) 

which was already considered in I. The second line of eq. (|25|) generates the already 
mentioned ambiguity: it is not uniquely defined since two different combinations of 
Green's functions can be used: 

61+2+3 (^16-26-3) = S 1+2+3 J l -^ e ^ + G { X^(E) (27) 



or 



5i+2+3 (6*16-26-3) = <5i+2-3 J l 4^e lEr]+ G { *\_ 2 j{E) (28) 

where the notations of the appendix (eq. (|5()|)) have been used. The second form 
gives identically zero in standard RPA and in r-RPA when covariance is forced in 
1^ (see appendix for its definition), whereas the first one gives a finite contribution. 
If eq. I)27[) is used, we obtain for the non covariant correlation energy: 

4o"r C) _ f 1 dp r _dP_ r idE iErj+ p 2 b 2 s 2 ( E + e pP \ 2I ( p 1) {E, P)Gp(E,P) 
V Jo P J (2Tr) d J (2vr) e 6 V 2e pP J l-^I p (E,P) 

(29) 

In 1 + 1 dimensions all the various contributions can be calculated using a Wick 
rotation (E 2 — P 2 — > —S) according to the method explained in I. 



E^l _ f 1 dp f°° dS 

~v~ ./„ 7./, 4^ 



E%t _ f 1 dp f°° dS p 2 b 2 s 2 I p (-S)G P (-S) 



V Jo p Jo ^ 6 1 _ glJ-S)) 
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Figure 1: Generalized mean field effective potential (left panel); generalized mean 
field effective potential plus the factorized term FT (right panel), for different values 
of the reduced coupling constant p, as a function of s. 
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4 Numerical results 

Before going to the discussion of the correlation energy, let us mention the following 
result. The generalized mean-field effective potential Eq(s) (vacuum energy as a 
function of the condensate s) presents a strong first order phase transition as shown 
on fig. ^ left panel (all numerical results are obtained with the reduced coupling 
constant p = 6/24/u 2 and fj, = 1). When we add the factorized term FT, this strong 
potential barrier disappears (fig. ^ right panel). The transition is practically, up to 
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Figure 2: Non covariant 3-body correlation energy £( 3nc ) calculated with eq. ((HEJ) 
(left panel: standard RPA, right panel: r-RPA), for different values of the reduced 
coupling constant p, as a function of s. 



an extremely small potential barrier, of second order nature with a critical coupling 
constant p c ~ 2.3. It is thus tempting to consider the sum Eq + FT as the true 
mean field energy for the r-RPA calculation because it contains all factorisable or 
reducible parts of the total energy. The fact that this term has a second order phase 
transition is very important in the following. The full energy is obtained by adding 
the pure interaction terms (H3+ : :} (diagrammatically, irreducible terms) which 
we will discuss below. 

We first consider the ambiguous term 

E (3nc) 

term using eq. (|33[) and we obtain 
the results shown on fig. [21 for standard and renormalized RPA. We can verify 
that this term is negligible at least in r-RPA. The maximum of this term is about 
respectively 10 and 100 times in s-RPA and r-RPA, smaller than the covariant 
contribution in the correlation energy shown on fig. El This demonstrates that the 
ambiguity linked to the non covariance problem is not so serious in r-RPA. Hence, 
in the following we take the reasonable option of putting this contribution to zero. 
The other term not considered in I, £K 4nc ) , is shown on fig. [31 in comparison with 
the covariant correlation energy £'( 3c ) + i?( 4c ) . We clearly see that E^ nc ^ which was 
ignored in I is sizeable especially for low values of the condensate s. The covariant 
effective potential (E cov = E Q + FT + E^ + E^) and the total effective potential 
(E tot = Eq + FT + E^ + E^ + £( 4nc )) in r-RPA are shown on fig. H For 
completeness we also show on fig. [21 the results obtained in standard RPA which 
were already given in I (in this calculation the RPA single particle propagator is 



0,2 - 



0,17b - 




0,1 0.2 O.J G.4 0.5 O.fi 0.7 0.8 

Condensate s 



o _p = _M - 

— - ———_ p _=_o 1 s_ 



——-—n=±.o 




-0.5 - 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Condensate s 



Figure 3: Left panel: non covariant 4-particle correlation energy E^ nc \ Right panel: 
covariant correlation energy 

£(3c) + £(4c) ( in r . RP A, for different values of p, as a 

function of s). 



replaced by the mean field one to avoid the divergence associated with instability). 
One sees that the RPA fluctuations in the standard RPA case are able to transform 
the strong first order transition of the Gaussian approximation into a second order 
phase transition. This comes from the fact that the attractive s-RPA correlations 
in the metastable region of the Gaussian approximation (p £ [0.2,0.5]) strongly 
reduce the potential barrier of the Gaussian effective potential. We see on fig. 03 
the evolution of the global minimum of the effective potential with increasing p. 
The fact that the transition is of second order is demonstrated by the fact that the 
condensate s is a continuous function of p and this is in agreement with the Simon- 
Griffiths theorem [S] which states that the phase transition in the Xcp 4 cannot be 
of the first order. However the critical coupling p c = 1.8 is different from both the 
lattice result [SIEB^J, Pc = 2.55, and the cluster expansion technique result [T^] . 
p c = 2.45. Although this result is not so bad, the standard RPA result cannot be 
really trusted since this method is spoiled by the instability problem. We obtained 
in I a preliminary result in r-RPA but keeping only the covariant pieces in the 
correlation energy. Here, the result is shown on the left panel of fig. 0] with the 
explicit incorporation of the FT term as explained in section 01 Although much 
less marked than in the Gaussian case, one again obtains a first order transition. 
There is nevertheless a slight progress with respect to the result obtained in I. The 
potential barrier is smaller and the critical coupling constant, p c = 1.9, is much 
closer to the lattice and cluster results. When the non covariant contribution E^ inc ' 
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Figure 4: Effective potential in r-RPA; Left panel: covariant potential E\J t = Eq + 
E^+E^+FT. Right panel: total potential E tot = Eq+eQ ) +£( 4c ) +FT. 



is added, the transition is of second order nature (see fig. HJ right panel). Indeed, 
this repulsive I?( 4nc ) decreases the correlation energy in the meta-stable region i.e., 
around the local maximum. Consequently the potential barrier is so weakened that 
it becomes negligible (about 10~ 3 ). Although the critical parameter is p c = 1.6, this 
constitutes an important result of this paper. The restoration of Simmon-Grifhth 
theorem indicates that r-RPA describes correctly the phase transition region. In 
addition the absolute value of the effective potential becomes very similar to the 
cluster effective potential. As an example for s = and p = 2, one gets the following 
results for the total energy. Ignoring the non covariant contribution one gets —0.14 
to be compared with —0.11 in the cluster calculation. When I?( 4nc ) is added, one gets 
—0.105. This is a very encouraging result and certainly further work is needed in 
direction of the fully self-consistent RPA. The problem of the covariance is certainly 
a key issue in that respect. One interesting possibility is to incorporate three-body 
excitations in the line of the work of ref. |13j . 

5 The linear— a model 

We now discuss another model, the linear-cr one. It can be seen as a Ginsburg- 
Landau effective Lagrangian for the SU(2) x SU{2) chiral symmetry (spontaneously 
broken in one direction). This model can give us physical insight in the chiral phase 
transition and may also describe far from the transition the dynamics of pionic 
systems. 
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Figure 5: Left panel: Covariant effective potential in s-RPA E\ c J t = E +E^+E^ c l 
Right panel: value of the condensate at the global minimum of E$. in s-RPA. The 
transition is of the second order, because this curve is continuous. 

The Lagrangian reads: 

C a = \d»od»o + - ^ (a 2 + n 2 ) -\{<r 2 + n 2 ) 2 + ca (35) 

where 7? is pseudo-scalar isovector field and a a scalar iso-scalar field, n corresponds 
to the physical pion and its chiral partner a may describe a mode associated with 
the amplitude fluctuation of the chiral condensate. 

This model can formally be seen as a generalization of the <J> 4 model for a N + 1 
dimensional multiplet (er, tt). It possesses an exact 0(N + 1) invariance if the pa- 
rameter c is zero. The ca piece of the Lagrangian describes the amount of explicit 
breaking of chiral symmetry in QCD. 

For the application of r-RPA we use notations similar to the 3> 4 model. We introduce 
s, the condensate in a field direction (the chiral symmetry being broken in one 
direction) and the fluctuating field a': a = a' + s with (a) = s (hence (a') = 0) and 
we omit the prime thereafter. We note Tl a and LT,,- the conjugate momenta of the 
fields a and n with usual commutation relations. 



The Hamiltonian reads: 



H(a, tt) = J dx 1 1 (Ul + (W) 2 + ( M g + Xs 2 )n 2 ) + (Va) 2 + ( M 2 + 3As V) 

+ As(j((T 2 + vr 2 ) + ^(a 2 + vf 2 ) 2 + a(/i 2 S + As 3 - c) + ^ 2 + ^ 4 - csj. 

(36) 

We define quasi-particle operators for quasi-pion and quasi-sigma: 



^ = y^ + yd^ and ^ = v^r^ + y2^ n ^ 

and we introduce the scalar densities jV Q = (7r a 7ia} and iV Q = (<j a <Ta). For what 
concerns the sigma operators, the indice a represents a momentum state whereas 
for the pionic operators it represents a momentum and isospin state. 

As in ref. we can introduce the r-RPA excitation operators in the pionic channel 
according to: 

— Y^ a _/3 b n> - a b a} 

This is equivalent, in the Green's function approach that we really use at variance 
with [2j, to calculate the RPA correction to the pion mass operator (E^) originating 
from the na RPA bubbles. 

With the same techniques used in the $ 4 theory, we derive the generalized mean 
field equations for pion and sigma modes. For N = 3 they read: 

4 = Mo + k l + A ^2(N a + 5M a ) + As 2 (37) 

a 

E 2 a = nl + kl + 3A ^{N a + M a ) + 3As 2 . (38) 

a 

The r-RPA inverse pion propagator is obtained as: 

G- l {u,k a ) =u? - el{uj) - Ej(w). (39) 
with the pion mass operator given by: 

E£(w) = 4AV 



,2 2 C(^) 



1-2A/J"H 

It contains an iteration of the ira bubble I ncT which is given by: 



or, in a form analogous to the RPA loop in $ model 



g/3 + -5/3' AZ/jg/j + NpiEpi _ ep — Ep> Mpep — NpiEp 
2e p Ep uj 2 -{ef3 + E pl Y~ 2e p Ep w 2 - (e p - E p ,)- 



T (41) 



Our preliminary goal was to show that the r-RPA fluctuations restore the Gold- 
stone theorem. With the use of the generalized mean field equations and using the 
expression of eq. I)39|) for u = 0, one can obtain the result: 



G- V, ka) =lo 2 - \i4 + As 2 + 3A ^(N a + M a )] - - ^(u = 0)) (42) 



(the last result (jl3"|) is obtained by using the gap equation dE/ds = which shows 
that the term in bracket in eq. (|42|) is just c/s). These expressions clearly show that 
the Goldstone theorem is satisfied in the chiral limit because the spurious mode 
uj = is allowed. 

As a preliminary conclusion, we underline this encouraging result: the r-RPA fluc- 
tuations can correctly treat the spontaneously broken symmetry and the spurious 
mode is obtained even if the covariance is lost. In particular it restores the Gold- 
stone theorem violated at the level of the mean field or Gaussian approximation [3] . 
This generalizes in r-RPA the result already obtained in the s-RPA formalism in [2] . 

6 Conclusion 

We have discussed in this article some problems encountered in the calculation of the 
effective potential in r-RPA and we have also presented some new numerical results 
for the A3> theory in 1+1 dimensions. We have shown that the ambiguity in the 
calculation of the 3-particle energy is only apparent in the sense that it is numerically 
very small. We have also shown that the incorporation of the so-called non covariant 
contributions in the effective potential significantly improves the description of the 
phase transition in the direction of lattice and cluster expansion results. The most 
important result of this work is that we found a way to take into account the different 
contributions of the RPA correlations that give us a second order phase transition. 

For what concerns the linear-cr model we have shown that the Goldstone theorem 
is explicitly satisfied despite of the covariance problem. A further work of interest is 
evidently to calculate the effective potential possibly at finite temperature to study 
the chiral phase transition. 




(43) 
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8 Appendix 



In this appendix we list the explicit expressions for the Green's functions. Here 
the momentum indices are represented by greek letters a, (3, 7. For this purpose we 
introduce various quantities: 

1 1 epMp + £p>Np> 



2ep 2ep, E - ep - ep> + irj 
ep - epi spNp - spiNpi 



Ispepi E 2 - (ep - Ep,) 2 + if} 
1 1 epMp + e'pNp< 



2ep 2ep> E + ep + ep< - ir) 



(44) 



We also introduce the loop integrals 

1 

V 



I a (E) = lg\E) + lg\E) + lg\E) . (45) 
In particular for a corresponding to the momentum P one has the explicit expression: 
I(E,P) = I a=p (E) 



$(d) (p-ki- k 2 



dk\ dk2 
(2vr) d 

£\ - £ 2 £iM - £2 A/2 



£\ + £2 £lM + £2 A/2 



2ei£ 2 E 2 - (£1 + £ 2 ) +irj 



2£i£ 2 E 2 - (£1 - £ 2 ) 2 + ir]_ 
For the one-particle Green's functions one obtains: 



(46) 



G aa't(E) = 8 a>a > 



E + £ a + 2e 



2£r. 



gaCg) 

— G a (E) 



—E + £ + Sa ( E ) 
G_ a] _ al {E) = 6 a , a , J 2 -^G a (E) 

G- a i a 't(E) = G a _ a >(E) = 5 a ^ a i 2^2 " G a (E), 



(47) 



where the full propagator is: 

G t (E) = 5 a ,a> G a (E) = 5 a , a , {E 2 - e 2 a - £ Q (£)) _1 . (48) 
The mass operator being given by: 

_ 11. 'f) . (49) 

For what concerns the 2p-lh and 2p-2p Green's functions we introduces indices i to 
label the destruction (creation) operators: 1 = (3, (3' '(ft, /?''), 2 = ((3, —(3^) sym ((ft , —13') 
and 3 = -ft, -ft(-(3, -ft). The results are: 

= w Sa -"-"' i-tuE) i^r) G ^ {E) 

G%>,-a(E) = <?%,(£) 

bs g I%,{E) f-E + e a 



^ Sa-B-B' — "I — G. , t (E) (50) 



b ^ Sa-P-P'Ipl'iE) (5 Q _ 7 _yJ^ / (£') 

Finally, we give the useful relations: 

E Ejf7 bS ( G f^W +G %',-*W) =^{E)G Mi {E) (52) 

i=l,2,3 PP' V 

and 

E E = 2 E + 6 E r^f 1 ^ 

ijel, 2,3/3/3' 77' a a 2 

2 s2 (I a (£)) 2 G (£) (53) 

These particular combinations are explicitly covariant if we forced it in the loop 
I(E,p)=I(E 2 -p 2 ). 
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